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Abstract 

There are non-radial null geodesics emanating from the shell focusing singu¬ 
larity formed at the symmetric center in a spherically symmetric dust collapse. 

In this article, assuming the self-similarity in the region filled with the dust 
fluid, we study these singular null geodesics in detail. We see the time evolu¬ 
tion of the angular diameter of the central naked singularity and show that 
it might be bounded above by the value corresponding to the circular null 
geodesic in the Schwarzschild spacetime. We also investigate the angular fre¬ 
quency of a physical held which propagates along the singular null geodesic 
and hnd that it depends on the impact parameter. Further, we comment on 
the non-uniformity of the topology of the central naked singularity. 
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I. INTRODUCTION 


Gravitational collapse is one of the most important issues in general relativity and in the 
theoretical research for this phenomenon, the Lemaitre-Tolman-Bondi(LTB) solution which 
describes the motion of spherically symmetric dust fluid has played an important role. This 
is a very simple analytic solution of Einstein equations but has substantial physical contents. 
The simplest but physically signihcant member of this solution is the case of a homogeneous 
dust sphere, which is called the Oppenheimer-Snyder solution. This solution is the hrst 
example to give an elegant picture of a black-hole formation and led to the so-called cosmic 
censorship conjecture; roughly speaking, this conjecture states that spacetime singularities 
except for the Big Bang are not visible for any observer if it is resulted from physically 
reasonable initial conditions [1]. However, this conjecture has not yet been proven although 
it plays very crucial roles in the proof of the singularity theorems [2-4] and the area theorem 
of black hole [5]. Rather LTB solution seems to be a candidate of the counterexample of 
this conjecture. Eardley and Smarr pointed out that a shell focusing singularity formed 
at the symmetric center in an inhomogeneous dust sphere can be naked, i.e., visible in 
principle [6]. After this discovery, several theoretical efforts revealed the genericity of the 
shell focusing naked singularity formation in the LTB solution [7-12]; there are radial null 
geodesics emanating from the central shell focusing singularity. 

The central naked singularity of the dust sphere seems to be a “point” since it forms 
from an inhnitesimal portion at the symmetric center of the dust sphere. However, it 
is not true. Recently, Mena and Nolan showed that there are non-radial null geodesics 
emanating from this central naked singularity [13], where non-radial means that it has a non¬ 
vanishing angular momentum. Deshingkar, Joshi and Dwivedi (DJD) studied numerically 
the trajectories of those geodesics within the dust fluid in detail [14]. More recently, Nolan 
and Mena studied various invariant quantities associated with this central naked singularity 
and further discussed the topology of this naked singularity [15]. 

It is the main purpose of this article to reveal how the central naked singularity is observed 
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for distant observers outside the dust sphere. For this purpose, we investigate both the radial 
and non-radial null geodesics emanating from the central naked singularity in the self-similar 
dust sphere. Assuming that the dust sphere has a hnite mass, we investigate the angular 
diameter of the naked singularity, or equivalently, the largest impact parameter of the null 
geodesics from the naked singularity; we show its time evolution. We also clarify the non¬ 
uniformity of topology of the central naked singularity for the self-similar case analyzed by 
Nolan and Mena [15]. 

We also study angular frequencies of physical fields which propagate along the null 
geodesics emanating from the central naked singularity. DJD showed that the redshift 
between the central naked singularity and any off-center observer is inhnite [14]. Hence 
even if physical helds with hnite frequencies are excited at the naked singularity, “physical” 
information of the naked singularity can not reach distant observers. However, it would not 
be natural that physical helds with finite frequencies are generated at the naked singular¬ 
ity where the energy density of the dust huid and spacetime curvature diverge. In reality, 
linear perturbation analyses of LTB solution show that the spacetime curvature associated 
to gravitational-wave perturbations indehnitely grows in the limit to the Cauchy horizon 
[16,17]. This means that the frequencies of the linear gravitational waves generated at the 
naked singularity will be inhnite even in purely classical dynamics. Physical helds excited 
by quantum ehects also have inhnite frequencies at the naked singularity [18-23]. Hence, 
in this article, we study the angular frequencies of the physical helds emanating from the 
central naked singularity with appropriately infinite frequencies. 

This paper is organized as follows. In Sec.H, we present the model of background space- 
time, which describes the gravitational collapse of a self-similar dust sphere with a hnite 
mass. In Sec.HI, we derive asymptotic behaviors of the null geodesics near the naked singu¬ 
larity and the basic equations for numerical calculations to see the intermediate behaviors of 
the null geodesics inside the dust sphere. In this section, we comment on the previous study 
about the topology of the naked singularity by Nolan and Mena. In Sec.IV, we explain the 
null geodesics in the exterior vacuum region. In Sec.V, we present the numerical results for 
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the apparent size of the naked singularity, i.e., the angular diameter of the central naked 
singularity, or equivalently, the largest impact parameter of the singular null geodesic. In 
Sec.VI, the frequencies of physical helds propagating along null geodesics emanating from 
the central naked singularity are investigated. Finally, Sec.VII is devoted for summary and 
discussion. 

In this article, we adopt the unit G = c = h = 1 and follow the sign conventions of the 
metric and Riemann tensors in Ref. [24]. 

II. COLLAPSING DUST SPHERE 

In this section, we show the model of background spacetime; the gravitational collapse 
of a self-similar dust sphere with a hnite mass. We assume that the dust fluid exists only 
within a hnite comoving region 0 < r < Tg, where is a positive constant. 

A. Interior Solution 

Dust huid elements move along timelike geodesics. We can adopt the tangent vector of 
these geodesics as a temporal coordinate basis. Then the line element is written in the form 

ds^ = —dt^ -\- ^ _l_ sin^ (2.1) 

1 + /(r) 

where a prime denotes the derivative with respect to the comoving radial coordinate r, and 
/(r) is an arbitrary function. In this coordinate system, the components of the 4-velocity 
of the dust huid elements are = (5f. Einstein equations and the equation of motion for 
the dust huid lead to 



(2.2) 

F' 


(2.3) 


where a dot means the temporal derivative, p(t, r) is the rest mass density of the dust huid, 
and F{r) is an arbitrary function. From Eq.(2.3), the arbitrary function F{r) is regarded as 


4 



the twice of a mass within the comoving radius r. Then Eq.(2.2) is regarded as an energy 
equation for a dust fluid element labeled by r and hence /(r) corresponds to the specihc 
energy of it. 

For simplicity, hereafter we assume the self-similarity for this spacetime; all the dimen¬ 
sionless variables are the functions of the self-similar variable x := t/r only. Then the 
arbitrary function F{r) is written in the form 

F{r) = Aor, (2.4) 


where Aq is a positive dimensionless constant. The specihc energy /(r) should vanish by 
this self-similar hypothesis. Then the solution of Eq.(2.2) is given in the form 


R = r 



(2.5) 


where using a freedom to re-scale the radial coordinate r, we have chosen an arbitrary 
function associated to the temporal integration so that R = r at t = 0. The dust huid 
element at r forms the shell focusing singularity at t = 2r/3^/Ao. 


B. Exterior Solution 


The exterior region is assumed to be vacuum and hence is described by the Schwarzschild 
geometry by Birkhoff’s theorem. Using the Schwarzschild static chart, the line element of 
the exterior region is written as 

J d2 

= —C{R)dT‘^ + + R? (dd^ + sin^ 9dip^\ , (2.6) 

l7(/x) ^ 

where the function C{R) is dehned by 


C{R) := 1 


2M __ Aors 

R R 


(2.7) 


The angular coordinates 9 and 99 are common to both the interior and exterior regions. 

The trajectory of the surface r = rg is a marginally bound timelike geodesic in the 
Schwarzschild spacetime by continuity. Thus the integration of the temporal component of 
geodesic equations and the normalization condition of the 4-velocity vector lead to 
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( 2 . 8 ) 


dT _ 1 

H ~ c{Ry 

dR _ 

Ht ~ “v 


(2.9) 


where t is the proper time which agrees with the interior temporal coordinate, and to derive 
the second equation, we have used the fact that the dust sphere is collapsing. 

From Eqs.(2.8) and (2.9), the trajectory of the surface is determined by 


dR 

dT 



( 2 . 10 ) 


The above equation is easily integrated and we obtain 


T-Td 


2 


(i? + 6M) — 2M In 


^/R — y2M 
\fR + y/2M 


where Td is an integration constant. 


( 2 . 11 ) 


III. NULL GEODESICS INSIDE THE DUST SPHERE 


By virtue of the spherical symmetry, it does not loose any generality to assume that the 
null geodesic is conhned within the 9 = ti/2 plane and hence we do so. Then the tangent 
vector of the future directed outgoing null geodesic in the LTB spacetime is expressed in 
the form 


k^ 


dt V 
dX R' 
dr _ - /2 

d\ ~ RR' 


dip I 

~d\^W 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


where A is the affine parameter and I is an integration constant which corresponds to the 
conserved angular momentum. Without loss of generality, we can assume / > 0. We call a 
null geodesic of / = 0 the radial one while a null geodesic of / > 0 is called the non-radial 
one. V should be larger than or equal to I so that /c* is positive and /c'’ is real. In this article. 


6 



we are interested in the null geodesics emanating from the central naked singularity; we call 
those the singular null geodesics. 

Instead of the affine parameter A, we adopt the comoving radius r to parameterize the 
null geodesic. Then the geodesic equations become 


dV 




dr 

dR 


1 


dr'^ 


^Ao(P2 - /2) + 

W- 


VV^ - P 



where we have introduced a new variable dehned by 


■ (W^ + 2r^d-«)/2^ _ 



and a is some positive constant which will be hxed in the following discussion. 


(3.5) 

(3.6) 


(3.7) 


A. Asymptotic Behavior near the Singularity 

As discussed by Joshi and Dwivedi [10], singular null geodesics should behave near the 
singularity as 

R —^IToV“ for r—^ 0, (3.8) 

where Wq is some positive constant. Assuming this behavior, we search for asymptotic 
solutions of both the radial and non-radial singular null geodesics in the neighborhood of 
the central naked singularity. 


1. Radial Singular Null Geodesics 

The asymptotic behavior of the radial singular null geodesic have already been studied 
well by several researchers [9-12]. However, we would like to again see it in detail since the 
knowledge about the radial singular null geodesics is useful for the analysis of non-radial 
singular null geodesics. 

For radial null geodesics, Eqs.(3.6) and (3.8) lead to 


7 



(3.9) 


(’^'0 - Aj''"’-''"”’''") (»o“ + 2r“<‘-«)/2) , 

r^o 3aWo^ A / \ / 

The above equation determines Wq and a. In order that the above equation has a positive 
root for hho, a should be 1/3 or unity. The solutions of a = 1/3 correspond to null geodesics 
from the origin before the formation of central singularity while the solutions of a = 1 are 
of our interest, i.e., the singular null geodesics. To see this fact, we rewrite Eq.(2.5) in the 
form 


t = 




1 - 


r J 


(3.10) 


Using the above equation and Eq.(3.8), we hnd that t —2 Wq/3AI^‘^ for r —0 along the 
null geodesic of a = 1/3. This means that the null geodesics of a = 1/3 start from the 
regular center r = 0 at t < 0 and hence these are not singular null geodesics. On the other 
hand, in the case of a = 1, f approaches to zero for r —0 along the null geodesic. Hence 
the null geodesics of a = 1 emanate from the central singularity r = 0 at t = 0. 

Choosing a = 1, Eq.(3.9) leads to 


yi/2 ^...3 21Uo(l-lUo3) 

Ao -H{Wo).- ^^3^2 • 


(3.11) 


1 /2 

The radial singular null geodesics have to satisfy the above equation. Since Aq is positive, 
0 < Wq < 1 should hold. We can easily hnd that there is a maximum of H(Wo) at 


/ 1/3 

Wo = := (-5 + 3^3) . 


(3.12) 


Hence in order that a singular null geodesic exists, Aq should be less than or equal to a 
critical value A^ dehned by 

Am := = {4 (26 - IsVs) . (3.13) 

We depict H{y) in Fig.l. In this article, we will focus on the case in which Aq is smaller 
than the critical value Am. Then form Fig.l, we can easily see that for a given value of Aq, 
there are two positive real roots of Eq.(3.11). We denote the smaller root by 1U_ and the 
other by Wc- 
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From Eq.(3.5), we find 


-l3{Wo)V for r—^ 0, 


(3.14) 


where 


mo) ■■= 


2 - lOlFo^ - IFo® 


31Fo3(lFo3 + 2) ■ 


(3.15) 


Since I3{Wq) 7 ^ 0 for 0 < Aq < A^, we obtain the asymptotic solution of V as 


V —Const. X r 


(3.16) 


From the asymptotic behaviors of W and V, it is better to introduce the following 
variables, 


P := 


AIF :=W- Wn. 


(3.17) 

(3.18) 


Since AW will vanish and P approaches to some hnite positive value Pq, for r —0, Eqs.(3.5) 
and (3.6) lead to 


dAW 

dlnr 

dP 


pAW, 

2Po(2 - 5W^) 


dlnr Wmo+2) 

Integrating the above equations, we obtain asymptotic solutions for r —0 as 


(3.19) 

(3.20) 


IF ^ IFo + IF,or^, 


^/i+ _ 6 ( 2 -5^0) _ 

r^\ Wo{2- lOWi - W^) j ’ 


(3.21) 

(3.22) 


where IF*o is an integration constant. 

Here again note that IF has to approach to IFq in the limit of r —> 0. This means that (3 
should be positive as long as IF*o does not vanish. The positivity of f3 leads to the following 
inequality 


0 < IFo < IFm. 


(3.23) 
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Hence for non-vanishing W^o, the smaller root Wq = W- of Eq.(3.11) has to be chosen (see 
Fig.l). 

The other root Wc of Eq.(3.11) is larger than and hence PiWc) is negative. For such 
singular null geodesics, hE*o has to vanish so that W approaches to Wc for r —0. This 
implies that W should be always equal to Wc by Eq.(3.6). This null geodesic is a generator 
of the Cauchy horizon. Substituting a = 1, / = 0 and W = Wc into Eq.(3.5), we can easily 
perform the integration and obtain the solution for V on the Cauchy horizon as 

V = Pcr-^% (3.24) 

where Pc is an integration constant and (3c ■= j3{Wc)- Since /3c < 0, "P vanishes in the limit 
r —> 0. 


2. Non-Radial Singular Null Geodesics ofV^oo for r ^ 0 


Let us consider non-radial singular null geodesics oiV-^ -|-cx) for r —0. The argument 
from Eqs.(3.9) to (3.23) is almost completely applicable to this case. However, because V 
has to diverge for r —0, /3 should be positive. Hence the smaller root Wq = W- of Eq.(3.11) 
must be chosen. 

The trajectory in the (r, </9)-plane is determined by the equation 

I dx i{w^ + 2y—^ 


dr SW^y/P^ - / V2/3-’ 

where (3- := j3{W-). The above equation becomes for r —> 0 as 

dif 1{WJ + 2)rl^--^ 

' 3,WJPo ■ 

Integrating the above equation, we hnd 

1{WJ + 2)r^- 




3/3_IE_3Po 


const. 


(3.25) 


(3.26) 


(3.27) 


Hence 99 has a hnite limit for r —0. As discussed by Nolan and Mena [15], this result 
means that the singularity from which the non-radial null geodesics of Wq = W- emanate 
is foliated by 2-sphere in accordance with Christodoulou’s argument [25]. 
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3. Non-Radial Singular Null Geodesics with Finite V\r=o 


Here, we investigate non-radial singular null geodesics with a finite limit of V for r —0. 
We denote this limit by Vq which should be larger than or equal to 1. However, we can easily 
see that there is no consistent solution of the behavior, V I for r —0 (see Appendix A). 
Hence hereafter, we assume Vq > 1. 

By the assumption, V is written in the form 

V = Vo + SV{r), (3.28) 

where 5V satisfies 


lim 5V = 0. 

r^O 


(3.29) 


Using Eqs.(3.8), (3.28) and (3.29), Eq.(3.6) becomes 

1 


Wn = lim 


3Q;Wn^ 


Wn 


^'■0 '0 




p 


rhP-o^) j I'pppS 2rid-")) . 


(3.30) 


In order that the above equation has a positive root Wq, a should be 1/3 or unity. By the 
same reason as in the previous cases, we find that a = 1 corresponds to the singular null 
geodesics, and hence we focus on it. Then the above equation leads to 




\ ' — 
— 


Vo (lUo' + 2) 

1 /2 

Since Aq' is positive, Wq is less than unity. 

From Eq.(3.5), we find 


(W 

dr 


A1/2 

'Wn^ 


V^o -v + 


Vo 


3IUn3 


(lUo^ + 2 


for r 


(3.31) 


(3.32) 


By the assumption oiV^Vo-, the right hand side of the above equation should vanish and 
hence following equation is obtained. 


,1/2 _VoWo{Wl + 2) 


(3.33) 
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From Eqs.(3.31) and (3.33), we obtain 


Pn = 


6/2(1 - Wo" 


(3.34) 


° 2-imo^ -Wo^' 

Since the right hand side of the above equation has to be positive, we obtain the same 
inequality as Eq.(3.23). It is worthy to notice that if Eq.(3.23) holds, is larger than 
Multiplying each side of Eqs.(3.31) and (3.33), we get 


Ao = l\Wo) := ^ (l - IFo=*) . (3.35) 

We also depict the function I(y) in Fig.l. From this hgure, we can easily see that there 
are two positive real roots of Eq.(3.35). Because of Eq.(3.23), we have to choose the smaller 
root which will be denoted by W+. Fig.l shows that W+ is always larger than W-. Hence 
the singular null geodesics of Wq = W+ are in the past of the singular null geodesics of 
Wq = W- in the sufficiently small neighborhood of the central naked singularity r = 0. 

By the assumption, we may rewrite W in the form 


W = W+ + SW{r), 

where 


(3.36) 


lim 5W = 0. 

r^O 


(3.37) 


Then in the limit of r 


0, Eqs.(3.5) and (3.6) are approximated by the linearized equations 


with respect to {SW, SV), 


where 


d 


d\nr 


Mn = 


Mi 2 — 


M21 — 


M 22 — 


^ 6W^ 

\6Vj 


2 - 10IF_, 


/ \ 

/ 

Mn Mi 2 


M21 M22 j 

V 

- IF +6 


6W 

5V 


3IF+3(IF+3 + 2) ’ 

^2(1 - IF+3)(2 - lOIF+3 - IF+6)3 
6^3 I Wl(W+^ + 2)2 ^ 

21 (2IF+3 + 1) 1^2(1 - IF+3) 
IF|^3(2 - lOIF+3 - IF+6) ’ 

2 - lOIF+3 - IF+6 

3IF+3(IF+3 + 2) 


(3.38) 

(3.39) 

(3.40) 

(3.41) 

(3.42) 
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Maiking the matrix Mij to be the diagonal form, Eq.(3.38) is written in the form, 

d 


dlnr 


f Si ') 17 0 W S, ') 

d -7y 


VS,; 




(3.43) 


where 


' Si ' 

1 ( 


2'yMi2 


iVfii + 7 Mi2 


\ ( AM/ \ 


y —Mil + 7 —Mi2 j 


SW 

ydVj 


(3.44) 


and 


_ J(4 - 8fE+3 - 5fE+6)(2 - lOfE+3 - W+^) 

7 :=/^ditM=- + 2) -' 


(3.45) 


Since Eq.(3.23) should hold, 7 is real. 
Solving Eq.(3.43), we obtain 


/sb / 


yS.; 


C'ir+^ 

C2r-^ 


\ / \ / 

where Ci and C 2 are integration constants. Since SW and 6V vanish for r 
be equal to zero. Hence we hnd that for r —0, 


W —^ fE+ + Wfor^, 
Mil - 7 


V 


V,- 




12 


fEtoM, 


(3.46) 
0 , C 2 should 

(3.47) 

(3.48) 


where W^o is an arbitrary constant. 

The asymptotic behavior of the trajectory in the (r, 99 )-plane is determined by 

d^ 1{W^ + 2) /(fE+3 + 2) 


dr SrWWV^ - P 3rW+^VW^ 
From the above equation, we obtain 

/(fE+3 + 2 ) 


for r 


0 . 




lnr +const. for r 


0 . 


(3.49) 


(3.50) 


3W+WVo^ - 

Hence we hnd that 99 —> —cx) in the limit r —0. As discussed by Nolan and Mena [15], 
this means that the singularity from which the singular null geodesics of IEq = 1 H+ emanate 
is topologically a “point” in the sense of Christodoulou [25]. As will be shown below, the 
singular null geodesics of IEq = 1 E+ have the largest impact parameter at each moment. 
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B. The Basic Equations for Numerical Calculations 


In order to see the intermediate behavior of singular null geodesics inside the dust sphere, 
we numerically solve Eqs.(3.5) and (3.6) with a = 1 by the 4th order Runge-Kutta method. 
By virtue of the existence of the homothetic Killing vector, we obtain an analytic expression 
for "P as a function of W (see Appendix B). However, we numerically integrate Eq.(3.5) and 
use the analytic expression for P to check the numerical accuracy. We will rewrite Eqs.(3.5) 
and (3.6) into appropriate forms for the numerical integrations below. 


1. Singular Null Geodesics ofV^oo for r ^ 0 


Here we consider singular null geodesics of P —cx) for r —0, or equivalently Wq = W_. 
From the asymptotic solutions (3.21) and (3.22), we introduce following variables. 


w := r^“, 

W^w) :=w-^AW. 


(3.51) 

(3.52) 


Then Eqs.(3.5) and (3.6) are written in the form, 

dP _2W, + {1 + 3p_) Pi + SAl^^wGi ^ 
dw~ 3f3-{Wo + wW,y 

dW, _ 2 (1 + 3/3_) W,P 2 + (2Wo + Ao^') P 3 + + 2 + wF2)G2 

~dw ~ 6/3_(Wo + wW*)3 

where 


(3.53) 

(3.54) 


Pi := + 6wW^W^ + WWoW^ + w^W^, 

(3.55) 

F 2 := 3W^W, + 3wWoW^ + 

(3.56) 

F3-.= 3WoW^ + wW!, 

(3.57) 

P 

C 

P (p + VP^ - ’ 

(3.58) 

P 

C ^ 

^ (p + VP2 - IV) VP2 - Pw^ 

(3.59) 
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Eqs.(3.53) and (3.54) can be easily integrated numerically from the origin tc = 0, or con¬ 
versely from the surface w = Wg := rj^- to tc = 0, since tc = 0 is not singular in these 
equations. 

It should be noted that W and P is not uniquely determined by its initial value Wq 
and Pq only. Eq.(3.54) shows that the boundary value of hh* is also necessary to determine 
them; we have denoted the value of IE*(0) by II4o- Since Wq should be equal to W-, the 
free parameters to specify the solutions for W and P are W^o, Po and the conserved angular 
momentum 1. However, here note that the solutions for W constitute a two-parameter 
family with respect to W^q and I/Pq. We can rewrite Eqs.(3.53) and (3.54) as the equations 
for P := P/Pq and hh*. By introducing I := I/Pq, Pq does not appear in these rewritten 
equations. The initial value of P is unity by dehnition. Hence the free parameters are lH*o 
and / in this new system of differential equations. 


2. Non-Radial Singular Null Geodesics with Finite V\r=o 


We consider the case of V —>£nite for r —0, or equivalently Wq = W+. In this case, the 
conserved angular momentum I is necessarily non-vanishing. From Eqs.(3.47) and (3.48), it 
is better to introduce following variables. 


V := r'’', 

W^{v) := v-^6W, 
V^{v) := v-^6V. 


Then Eqs.(3.5) and (3.6) are written in the form, 

dPt _ 1 


dv 

dW^ _ 1 

dv 6'yWH‘^ 


{W{W^ + 2) - 3Ay VP2 - /2 - 3-fW%v} , 
2IE(1 - W^) - - 6-fW^W^v } . 


Eqs.(3.47) and (3.48) impose boundary conditions on and as 


(3.60) 

(3.61) 

(3.62) 


(3.63) 

(3.64) 
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(3.65) 

(3.66) 


l^t(O) = W^to, 

IVl 12 

where hhj-o is an arbitrary constant. The right hand sides of Eqs.(3.63) and (3.64) seem to 
be singular at n = 0, but it is not true. Taking into account Eqs.(3.65) and (3.66), we can 
verify that the right hand sides have hnite limits for n —0. 

Once the initial value Wfo is given, the initial value of is uniquely hxed by Eq.(3.66). 
Hence the solutions of W^{v) constitute an one-parameter family with respect to Wfo (the 
solutions of 1T| seem to constitute a two-parameter family with respect to Wfo and /, but it 
is not true). In order to see this fact, let us introduce following variables 


^o:= 



n 

I 

h. 

I ’ 


6(1 - IE+3) 

^ 2-10IE+3-IE+6’ 


(3.67) 

(3.68) 


and rewrite Eqs.(3.63) and (3.64) into the equations for these variables. Then the parameter I 
disappears in the new equations. This means that the solutions for and Hj do not depend 
on 1. Since Vo is not an arbitrary constant by Eq.(3.67), the remaining free parameter is the 
initial value W-^o only. 


IV. NULL GEODESICS IN THE EXTERIOR REGION 


The components of the null geodesic tangent on the equatorial plane 0 = 7r/2 in the 
exterior Schwarzschild static chart are given by 


k^ 

k^ 


dT 

lx 

dR 

lx 


00 


C{R)' 


= 00 


\ 


dd ^ 

lx ~ ’ 


C{R) 


dip I 

IX^R^' 



(4.1) 

(4.2) 

(4.3) 

(4.4) 
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where u and I are integration constants which correspond to the angular frequency measured 
by the static observer at the future null inhnity and to the conserved angular momentum, 
respectively. 

In order to obtain a relation between the components of the interior and of the exterior 
regions, we introduce a tetrad basis whose components with respect to the interior LTB 
solution are given by 


= (»■ O' O) • 

= (•>■ O' O' ' («) 


Note that e(q^ agrees with the 4-velocity of the dust fluid element and e(r)^ is normal to 
e(q^ in (t,r)-space. At the surface of the dust sphere r = r^, the components of e(p^ with 
respect to the exterior Schwarzschild static chart are obtained by Eqs.(2.8) and (2.9). The 
components of e(r)^ with respect to the exterior static chart are obtained as an outward 
normal to e(t)^ in (T, i?)-space. Hence the components of the tetrad basis with respect to 
the exterior static chart are given by 



Using the above expression, the continuities of e(t)^k^ and e{^r)^k^ lead to the relation between 
the components in the exterior and in the interior coordinate systems at the surface of the 
dust sphere; 




1 



(4.13) 
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(4.14) 





+ B!k\ 


Therefore, we obtain the angular frequency for the observer at the null infinity as 



(4.15) 


The (/^-components k‘^ in exterior and interior coordinate systems agree with each other and 
hence I is common. 

We consider the situation in which a distant static observer at some given areal ra¬ 
dius R = Ro detects the null geodesics. There is the invariance of the time translation, 
T —0 T-|-const., in the exterior region. Using this symmetry, we set the origin of the time 
coordinate T so that the observer intersects the Cauchy horizon at T = 0. The future 
directed outgoing radial null geodesic in the exterior Schwarzschild region is determined by 


dT _ 1 

dR ~ C{R) 

We can easily integrate the above equation and obtain 


(4.16) 


T-Tng = R + 2M\n\R-2M\, 


(4.17) 


where T^g is an integration constant. Since the Cauchy horizon is generated by radial null 
geodesics which go through R = R^ at T = 0, the trajectory of the Cauchy horizon should 
be written in the form 


T — Tdfj(.R) •— R — Rq T 2A4 In 


R-2M 


(4.18) 


- 2M 

The Cauchy horizon intersects the surface of the dust sphere when its areal radius is given 
by 

2MWr^ 


R = Rr^ := r.Wc = 


An 


(4.19) 


Hence the exterior time coordinate at the intersection R = R^s is given by 


T — Tcs Tch(.Rcs)- 


(4.20) 
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The trajectory of the surface of the dust sphere is given by Eq.(2.11). The integration 
constant in this equation is determined by the condition in which the areal radius of the 
dust surface is equal to i?cs at T = Teg. Then we obtain 


Td 



{Res + 6M) + 2M In 


— V2M 
V Res + \/2M 


(4.21) 


To obtain non-radial singular null geodesics in the exterior region, we integrate the null 
condition 


^ __ R _ 

dR C{R)^R^ -{I/uyC{R)' 

between the surface of the dust sphere R= Rg and the observer R = Rq. 


(4.22) 


V. APPARENT SIZE OF THE NAKED SINGULARITY 


First we consider the angle 5 between the radial direction and the singular null geodesic 
at the observer R = R^ (see Fig.4). This angle 5 is dehned as follows. We introduce an 
orthonormal basis whose components with respect to the Schwarzschild static chart are given 
by 


(e(T)'")sch — ( ! -, 0, 0, 0 I , 

(5.1) 

VV^(^o) ) 


(e(ij)^)sch = ^0, yC(Ro), 0, 0^ , 

(5.2) 

(e((9)^)sch = ^0, 0, —, 0^ , 

(5.3) 

(e^ )sch- (o, 0, 0, 

(5.4) 


Then the angle 5 is given by 


tan 5 = 




I 


C{Re 


e(^Rfk^ u;\\Re^-{l/u;yC{Ro)' 

We can easily see that for sufficiently large Ro, the angle S becomes as 


(5.5) 




(5.6) 
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where 6 is a constant which is regarded as the impact parameter of the singular null geodesic 
with respect to the central naked singularity, 


b := lim Ro6 = 

R^^oo ijj ' 


(5.7) 


The angle S depends on the location of the observer. On the other hand, the impact 
parameter b does not and hence it is more convenient to characterize the “size” of the naked 
singularity. 

Eq.(4.2) is regarded as an energy equation of a test particle with a specihc energy 1/2, 

2 


1 (dRy ,, 1 

2 U = 


(5.8) 


where \ = u\ and 


V{R- b) 


b^C{R) 

2R^ 


(5.9) 


The maximal value of the effective potential V is 6^/54M^ at i? = 3M. Hence if the impact 
parameter b is smaller than 3v^M, there is no forbidden region for the motion of the test 
particle and hence the null geodesic can escape to inhnity if it is initially outgoing. On the 
other hand, in the case of & > 3\/3M, it is non-trivial whether the singular null geodesics 
can go away to the inhnity. When b is equal to dv^STf, the null geodesic of the circular orbit 
R = 3M is possible although such an orbit is unstable. 

Numerically integrating geodesic equations derived in Sec.Ill, we can draw trajectories 
in the {w, IT)-plane. As an example, we depict the radial singular null geodesics of IT_ = 
Y^d/lO in Fig.2. By Eq.(3.11), Aq is about 0.288. The radial singular null geodesic of 
ITo = W- approaches to the Cauchy horizon hh = 144 in the limit of H4o —^ oo. On the 
other hand, in Fig.3, we depict the trajectories of the non-radial singular null geodesics of 
Wq = W±. As in Fig.2, we choose IF_ = ^3/10. In this case, IF+ is about 0.361. All 
the null geodesics in this hgure are detected by the observer at Ro = lOOM at the moment 
T = 5M. Hence these null geodesics constitute an one-parameter family with respect to the 
impact parameter b. The lower most curve corresponds to the radial singular null geodesic 


20 



6 = 0. There is a critical value 6crit for the impact parameter b (6crit is about 0.460 x in 
the case of this hgure). The null geodesics of 6 > 6crit do not start from tc = 0 and hence are 
not singular null geodesics; W of these null geodesics becomes larger than ITc near tc = 0. 
On the other hand, the null geodesics of 6 < 6crit start from the singularity tc = 0. For 
b < 6crit, W is equal to W- at tc = 0, while in the case of 6 = 6crit, W becomes W+ at 
tc = 0. Hence the singular null geodesic of ITo = W+ has the largest impact parameter(see 
also Fig.5). 

The angular diameter A of the central naked singularity for an observer at i? = i?o is 
the twice of the angle 6 of the singular null geodesic with the largest impact parameter 6crit, 


A = 2 arctan 


-'crit 


C{Ro 


(5.10) 


In order to study the temporal behavior of the largest impact parameter b = 6crit, we 
integrate Eqs.(3.63) and (3.64) from u = 0 to u = rs''' for various boundary values W^q. We 
set I to be unity. After the integration to u = rg'*' is completed, we obtain u by Eq.(4.15); we 
denote u of the null geodesic of 6 = 6crit by cUcrit- Then we integrate Eq.(4.22) until R = Ro 
and as a result, we obtain Tq. Here note that Tq is a function of the initial value W^q. We 
search for the initial value so that some given To is realized as a result of the numerical 
integration. Then we obtain the largest impact parameter by 6crit = ^ait for given Tq. 

In Fig.6, we depict the largest impact parameter 6crit as a function of the exterior time 
coordinate To at the observer R^ = lOOM. Since Ro is sufficiently larger than 2M, the 
angular diameter A of the central naked singularity is almost equal to 26(.rit /Ro. 

From Fig.6, we see that the largest impact parameter 6crit approaches to 6 b := Sy/SM 
asymptotically. Hence our numerical results strongly suggest that the largest impact pa¬ 
rameter 6crit is bounded above by 6 b. The angular diameter A is also bounded above by 


Ab = 2 arctan 


6b 


C{Ro) 

Ro^ - 6b2C'(To) 


(5.11) 


First, at the Cauchy horizon To = 0, the central naked singularity is observed as a point, i.e., 
forit = 0, or equivalently A = 0. Then gradually the largest impact parameter 6crit grows 
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and approaches to 6 b; correspondingly, the angular diameter approaches to Ab- 


VI. FREQUENCIES OF PHYSICAL FIELDS FROM THE NAKED 

SINGULARITY 


The redshift zq between the central naked singularity and a timelike observer with a 
4-velocity is given by 


Zq + 1 = - -lim ki.u''. (6.1) 

r^O 

Since ki,u‘' = —V/R —oo for r —0 for hnite kfj,v^, the redshift zo is inhnite. From 
this fact, it seems to be likely that physical information of the naked singularity cannot 
arrive at any observer. However it may not be necessarily true. Since the energy density 
and spacetime curvature are inhnite at the naked singularity, physical helds generated there 
might also have inhnite frequencies. Hence even if the physical helds suher the inhnite 
redshift, it is a non-trivial issue whether their frequencies observed at inhnity vanish or not. 

We do not exactly know what and how amount of physical helds are generated in the 
neighborhood of the central naked singularity where quantum ehects of gravity might be 
important. However we may say that the physical held which can arrive at the distant 
observer from the dense region around the naked singularity might be graviton due to the 
weakness of its coupling to other helds. Further, it might not be so terrible to assume that 
the frequency of the graviton is order of ^/p which is the inverse of the free fall time. From 
Eq.(2.3), along a singular null geodesic, the density p behaves as 


3WoHI - 1 

27r(H'i“ + 2)3 



( 6 . 2 ) 


Hence hereafter we will focus on a physical held which is generated at R = L with a comoving 
angular frequency ul = l/L. Of course, L is assumed to be sufficiently smaller than M. It 
should be noted that the angular frequency of the physical held propagating along such a 
null geodesic becomes inhnite at the central naked singularity. 
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A. On the Cauchy Horizon 


First, let us consider the angular frequency of the physical held which propagates on the 
Cauchy horizon. As mentioned, the Cauchy horizon is R = Wc‘^r. Hereafter we will denote 
the comoving radius of the singular null geodesic at i? = L by ri- Then is equal to L/Wc^ 
on the Cauchy horizon. From Eq.(3.24), the angular frequency at i? = L is given by 


~ ■ 

As mentioned, we assume that the angular frequency ojl is 1/L aX R 
the above equation, we obtain 


Pc = 


L 


0c 


(6.3) 

L and hence from 


(6.4) 


From Eq.(4.15), the angular frequency for a static observer at inhnity is given by 


_ 1 / L 2fF,V2(i _ 

2M\2m) I (fFe + 2)2 J \ ^ Wc + 2 
It should be noted that since /3c is negative, a; diverges in the limit of L —0. 
We hnd that in the limit of HA —1, the angular frequency u behaves as 



u 


1 

L 


(6.5) 


( 6 . 6 ) 


This is because HA = 1 corresponds to no dust-sphere case Aq = 0. Hence in the case of 
sufficiently small Aq, the angular frequency u observed at the inhnity is almost equal to the 
initial value ul at R = L. However it should be noted that small Aq does not necessarily 
corresponds to small total mass M. Because of M = Aorg, if the comoving radius of the 
dust surface is very large, M can also be large. On the other hand, in the limit of HA —^ HAi, 
the angular frequency becomes as 

0^ ^ ll +(~ 0.243M-A (6.7) 

2MWra \ V HAi / J ^ 

The angular frequency u does not depend on L in this limit. 
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B. In the Future of the Cauchy Horizon 


Now let us consider singular null geodesics in the future of the Cauchy horizon. The 
angular frequency of the physical held propagating along the singular null geodesic of the 
largest impact parameter b = 6crit is obtained by integrating Eqs.(3.63) and (3.64) from 
n = 0 to n = Ts'^. For the exterior region, we integrate Eq.(4.22) until R becomes some given 
Ro to obtain T at R = Ro- Then the initial value of is determined iteratively so that T 
becomes some given Tq at R = Ro- The initial value of V is unity, i.e., Vo = 1, in order that 
ul = l/L. This condition and Eq.(3.34) lead to 


/ 


^crit • 


2 - lOfE+3 - 

\l 6(1 - fE+3) 


( 6 . 8 ) 


As a result of the numerical integration, we obtain the angular frequency cUcrit and the impact 
parameter 6crit = Icnt/^cnt- Here note that the angular frequency cUcrit does not depend on 
the assumed value of L. 

In order to study the angular frequencies for the distant static observer in the case of 
the null geodesics of Wq = W-, we had better integrating the geodesic equations from the 
observer toward the central naked singularity. First, we integrate Eq.(4.22) from the observer 
R = Ro toward the dust sphere R = Rg. The initial condition is the given To- We set u 
to be unity and hence the impact parameter is equal to / which should be smaller than bent 
obtained in the above. After the null geodesic reaches R = Rg, we integrate Eqs.(3.53) and 
(3.54) toward the central naked singularity tc = 0. The initial condition for these equations 
are obtained by Eqs.(4.13) and (4.14) as 


P{Wg) 

and 


Once we obtain P at R 
observer at R = Ro as 


WgRg 


W.{wg) = 



2M \ 




1 RqRs 

Wg\\ 

2M 


PC{Rg) 

2 


(6.9) 


W- 


( 6 . 10 ) 


: L, we obtain the observed angular frequency ojo for the static 
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( 6 . 11 ) 


- 

P{wl)^C{rS 

where Wl = r. 

In Figs.7 and 8, we depict the angular frequency Uo for the static observer at Ro = lOOM 
as a function of the impact parameter b. Fig.7 corresponds to the case of L = 10“^M while 
Fig.8 is for L = 10“^M. Since Ro is sufficiently large, the angle S is almost equal to b/Ro- 
The background model is chosen so that IF_ = y^S/lO. There are four curves; these are 
the data at Tq = 5M, lOM, 15M and 20M, respectively. The largest impact parameter 
&crit normalized by 6 b for each moment is about 0.460, 0.736, 0.882 and 0.951, respectively. 
The angular frequency is monotonically increasing with respect to the impact parameter b 
or equivalently, the angle S, while it monotonically decreases with the lapse of time. The 
angular frequency at 6 = 0 will vanish in the limit of Tq —> cxd, since in this limit, the radial 
null geodesic agrees with the generator of the event horizon whose redshift is inhnite. 

The angular frequency cUcrit is depicted as a function of the time To in Fig.9. This hgure 
suggests that ojait approach to some hnite value in the limit of T oo. Hence the naked 
singularity will be observed as a ring for sufficiently late time. 

VII. SUMMARY AND DISCUSSION 

We studied non-radial null geodesics emanating from the central naked singularity formed 
by the gravitational collapse of a self-similar dust sphere. The existence of the non-radial 
singular null geodesics means that the central naked singularity is observed as a disk if it 
emits null geodesics of all possible impact parameters. The observed singular null geodesic 
with the largest impact parameter b = bent gives the angular diameter A of the naked singu¬ 
larity. The angular diameter A is time dependent; it monotonically grows and approaches 
to the value 3t/3M/Ro for the observer at i? = To 3> M. The asymptotic value of the 
angular diameter comes from the geometry of the exterior Schwarzschild region. 

As shown by Nolan and Mena [15], in the case of the self-similar LTB solution, the 
topology of the central naked singularity is not uniform, being point-like on one region 
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and spherical on another. In this article, we have shown that the point-like region of the 
central naked singularity is a portion from which singular null geodesics of the largest impact 
parameter b = bent emanates. In other words, the outer boundary of the observed disk 
constitutes the null geodesics from the point-like portion of the central naked singularity. 
On the other hand, all the other singular null geodesics emanate from the spherical region 
of the central naked singularity. 

We also studied the angular frequencies of the physical helds which propagate along 
the singular null geodesics. It is well known fact that the redshift between this naked 
singularity and any timelike observer is inhnite. Hence even if some physical held with 
a hnite frequency is excited at the central naked singularity, it cannot carry its energy 
from the central naked singularity to the outside. From this result, it seems to be likely 
that no physical information comes form the central naked singularity to any observer. 
However, it might be terrible to assume that some physical held with a finite frequency 
are generated at the naked singularity at which the energy density of the dust huid and 
spacetime curvature are infinite. Graviton might be only one physical held which can 
carry the physical information to distant observers from the central naked singularity since 
it might not be scattered nor absorbed even in the situation with the extremely high energy 
density like as in the neighborhood of the central naked singularity. Its frequency might be 
order of ^/p, that is, the inverse of the free fall time. Since the energy density p near the 
central naked singularity is proportional iol /along the singular null geodesic, the angular 
frequency of this held generated aX R = L are order of 1/L. Hence in this article, we have 
considered the null geodesics which start at R = L with the comoving angular frequency 
ul = l/L and then have investigated the angular frequency Uo observed at i? = i?o 3> M. 
On the ground of dimensional analysis, if the helds are excited by the quantum gravitational 
ehects, L will be order of the Planck length. 

As a result, we found that the angular frequency is a function of the impact parameter 
b, or equivalently the angle S between the radial direction and the null geodesic. At the 
center <5 = 0, the angular frequency u is the smallest, while it is the largest at the boundary 
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6 = A/2. Hence, if such null geodesics are really emitted, the naked singularity is observed 
as a disk-like rainbow. 

In the limit of T —> cx), the angular frequency a; at 5 = 0 vanishes while it does not for 
0 < S < A/2. The singular null geodesic at 5 = 0 in the limit of T —cx) is the generator 
of the event horizon and hence suffers the inhnite gravitational redshift. By virtue of this 
gravitational redshift of the central naked singularity, we might observe the gravitons in a 
wide range of the angular frequency. Hence all the gravitational wave detectors (TAMA 
[26], LIGO [27], VIRGO [28], GEO [29], LISA [30] and DEGIGO [31]) now in operation or 
in plan are able to detect the graviton from this naked singularity if sufficient amount of 
gravitons are excited there. 

There are several researches for more general spherically symmetric systems and those 
revealed that non-vanishing pressure does not necessarily prevent the formation of the central 
shell focusing naked singularity [32-35]. In these example, there might be non-radial null 
geodesics emanating from the central naked singularity. The investigation of these is future 
work. 

Acknowledgements 

We are grateful to T. Harada and the colleagues in Department of Physics, Osaka Gity 
University for useful discussion. KN also thanks to the workshop of DEGIGO held in 
National Astronomical Observatory of Japan in May, 2002. 

APPENDIX A 

In this appendix, we show that there is no singular null geodesic of V ^ I for r —0 in 
the self-similar LTB spacetime. 

Proof. First we assume P —*> / in the limit r —0. Hence, V is written in the form 

V = l + SV{r), (Al) 

where 6V satishes 
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lim 5V = 0. 


(A2) 


Then in the limit r —0, Eq.(3.6) becomes 


Wo^ = lim ■ 


lEn- 


llAor 


1—a 


25V 




(A3) 


»’^o 3 q;1To^ 

Since the right hand side of the above equation should be hnite, a is less than unity. There 
fore, in the limit of r —0, Eq.(3.5) leads to 


d5V 

dlnr 

Integrating the above equation, we obtain 


3‘ 


(A4) 


5V —^ - Inr for r —0. 

3 


(AS) 


This behavior is in contradiction to the assumption (A2). □ 


APPENDIX B 

In the self-similar LTB solution, there is a homothetic Killing vector which satisfies 

£^9cLb 9ab' (® 1 ) 

In this spacetime, its components are given by 

e = (t,r,0,0). (B2) 

We can easily see that is constant along a null geodesic with tangent k^. By 

virtue of this fact, we obtain an analytic expression for V [9]. 

The self-similar variable x is expressed by W as 

"= 3 A^(' 

Hence R' can be also expressed by W as 

IK3 + 2 


R' = W^+ r{drW\ 


3IE 


(B4) 



The conserved quantity associated with the homothetic Killing vector is 


= tk‘ - rR’\r 


(B5) 


From the null condition = 0 and Eq.(B5), we obtain for non-vanishing C^, 

C ( /-1 

V = ± flVl - {x‘ - R’^)PHCiWY\^ (B6) 

When = 0 but I does not vanish, the positive root of V is 


-x^ 

For the case of = 0 = /, the solution of Eq.(3.54) is 
can be easily integrated as 


(B7) 


0. Using this solution, Eq.(3.53) 


V=—. (B8) 

w 

We numerically integrate Eqs.(3.53) and (3.63). Then the solutions (B6), (B7) and (B8) are 
used to check the accuracy of the numerical integrations. In our numerical calculations, the 
relative error estimated by using these analytic solutions for V is less than 10“^. 
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FIGURES 



FIG. 1. The functions H{y) and I{y) are depicted. The roots of Eqs.(3.9) and (3.11) are also 
depicted for a given Aq which is only one free parameter of the self-similar dust sphere. 



FIG. 2. The trajectories of radial singular null geodesics I = 0 are depicted in {w, IF)-plane. 
From this figure, we can see that in the limit of IF*o ^ oo, the singular null geodesic approaches 
to the Cauchy horizon W = Wc- We also plot the apparent horizon F/R = 1, or equivalently 
W = ill this plane. 
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w/M^ 

FIG. 3. The trajectories of singular null geodesics with various impact parameters b are depicted 
in (w, IF)-plane. These null geodesics intersect with each other at (T, R) = {5M, lOOM) in the 
exterior region. The lower most curve corresponds to the radial null geodesic / = 0 (6 = 0). The 
dashed curve is the non-radial singular null geodesic with the largest impact parameter b = 6crit- 
The solid curves below this dashed curve are all the singular null geodesics while those above this 
curve are not. 


Null Geodesic 



FIG. 4. The schematic diagram of the trajectory of a non-radial singular null geodesic. 
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FIG. 5. We depict the conformal diagram of the spacetime with a self-similar dust sphere. The 
singular null geodesic of 6 = 6crit emanates from the beginning of Cauchy horizon in this diagram. 
The observed null geodesics of 6 > 6crit do not go through the origin r = 0. 
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FIG. 6. We depict the impact parameter ftcrit as a function of time To at the observer 
Ro = lOOM. The crosses are the data for the background self-similar dust chosen so that 
W- = a/2/10 while the squares correspond to the data of the background with W- = a/3/10. 
The impact parameter ftcrit grows monotonically and approaches to the value 6 b = 3\/3M. 




6/6b 

FIG. 7. The angular frequencies ujo for the observer R = Ro at T = 5M, lOM, 15M and 20M 


are depicted as a function of the impact parameter 6. The null geodesics are assumed to start at 
R = W~^M with the comoving angular frequency ujl = The background self-similar dust 


is chosen so that W- = a/3/10. The angular frequency lVq is a monotonically increasing function 


of 6. 
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FIG. 8. The same as Fig.7 but the null geodesics are assumed to start at R = 10~^M with the 
comoving angular frequency lol = 
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FIG. 9. We depict the observed angular frequency Wcrit of the null geodesic with the largest 
impact parameter b = 6crit as a function of time To at Rq = lOOM, where ujl is set to be 
The crosses are the data for the the background self-similar dust chosen so that W- = a/2/10, 
while the squares correspond to the data of the background with IT_ = a/3/IO. The observed 
angular frequency Wcrit approaches to some finite value for T ^ oo. 
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